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Abstract 

0^ In the context of Brans-Dicke theories, eternal inflation is described in such a way 

1^ that the evolution of the inflaton field is determined by the value of the Planck mass 

^ in different regions of the universe. The Planck mass is given by the values of the 

Brans-Dicke field, which is coupled to the scalar curvature in the Lagrangian. We first 
calculate the joint probability distributions of the infiaton and Brans-Dicke fields, in 
order to compute the 3-volume ratios of homogeneous regions with arbitrary values 
^ of the fields still undergoing infiation with respect to thermalized regions. From these 

O volume ratios one is able to extract information on the values of the fields measured 

by a typical observer for a given potential and, in particular, the typical value of 
a the Planck mass at the end of infiation. In this paper, we investigate volume ratios 

using a regularization procedure suggested by Vilenkin, and the results are applied 
to powerlaw and double- well potentials. The spectrum of density fiuctuations is 
calculated for generic potentials, and we discuss the likelihood of various scenarios 
that could tell us whether our region of the universe is typical or untypical depending 
on very general bounds on the evolution of the Brans-Dicke field. 
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1 Introduction 



It is widely accepted within the standard framework of quantum cosmology that constants of 
nature take different values in different regions of the universe (Linde f 986,1 989,1 990a; Cole- 
man 1988; Vilenkin f995a,b; Garcia-Bellido & Linde 1995). One of the chief goals of quantum 
cosmology is to predict the probability distributions for these constants (e.g. the gravitational 
constant G, the cosmological constant, the density parameter 0, etc) from the wave function of 
the universe. By predicting the values of the constants of nature for a "typical" region of the 
universe from first principles, one would in principle have a means to verify whether we do indeed 
inhabit a typical region of the universe or not by comparison with the observational data. This 
is however a difficult undertaking, and it is customary to adopt the reverse approach, i.e. that of 
assuming that it is likely that we do indeed inhabit a typical region of the universe and, therefore, 
our predictions for the values of the constants of nature should be consistent with the data in 
the observable universe. This is a form of "principle of mediocrity" that has been enunciated in 
terms of anthropic arguments in the literature (Carter 1983; Barrow & Tipler 1986; Rees 1993; 
Albrecht 1995). 

In models of eternal infiation, it is non-trivial to derive the probability distributions for the 
fields directly from the wave function of the universe (Linde 1995; Vilenkin 1995b), which is 
based on the Euclidean approximation to the sum over histories. This approach is not as yet 
entirely understood and yields unrenormalizable probability distributions for the fields that are 
diflicult to interpret (Hartle & Hawking 1983; Vilenkin 1984; Hawking 1987; Linde 1984,1995). 
Alternatively, the Fokker-Planck or diffusion equation enables us to compute these probability 
distributions based on the idea that the infiaton field evolves stochastically and thus takes dif- 
ferent values within a given region of the universe in a chaotic manner (see e.g. Starobinsky 
1984,1986; Linde 1986). The evolution of the field is governed by Langevin-like equations of 
motion, where the classical motion of the field is perturbed by its Brownian motion in the form of 
quantum jumps around the classical slow-roll solution. As a result, regions where the field rolls 
down the potential below the end-of-infiation boundary are thermalized, whereas those in which 
infiation still takes place are subdivided in further regions where the field takes different values, 
some of which thermalize and others continue infiating (Linde 1986,1987; Goncharov et al. 1987; 
Garcia-Bellido 1994). In this scenario, provided one starts out from sufficiently homogeneous 
initial conditions, at any given time there will be a fraction of the total volume of the universe 
that has thermalized, V*, and the rest still undergoes infiation. Chaotic infiation predicts that 
this process is eternal and it can be implemented for a wide class of potentials (see e.g. Vilenkin 
1983; Linde 1986). Also, it is found (Linde & Mezhlumian 1993) that the fraction of the physical 
volume where the infiaton field has a certain value a attains a constant value in the t — )• oo limit. 
This is what we will refer to hereafter as the concept of "global stationarity" . The implications 
of global stationarity for structure formation are strong as we shall see in the following. 

In order to compute the probability distributions of the constants of nature or the spectrum 
of density fiuctuations in different regions of the universe, it is necessary to compare the physical 
volumes of the hypersurfaces on which the fields take certain constant values. Provided that a 
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sufficient amount of time has ellapsed since the onset of inflation, so that global stationarity can 
be assumed, the fraction of the physical volume of the universe where a constant of nature has 
an arbitrary value will be constant. In the case of the gravitational constant G for example, we 
can investigate the spectrum of fluctuations in thermalized regions of the universe where G takes 
different values. The value of G is related to the Planck mass, and we expect this to be slowly 
varying in different regions of the universe, as described by a Brans-Dicke (BD) theory of gravity 
(see e.g. Brans & Dicke I96I; Brans 1962; Dicke 1962). Regions of the universe where the value 
of G is greater will experience a quicker process of galaxy formation and cosmic structures will be 
clumpier for the same value of the density parameter 0. Inflation is brought in this scenario by 
models of extended inflation (see e.g. La & Steinhardt 1989; Linde 1990b) as well as inflationary 
scenarios that incorporate in the Lagrangian a BD gravity theory (Linde 1990a, 1994a; Laycock 
& Liddle 1994). 

In chaotic inflation with BD gravity (Linde 1990b; Garcia-Bellido 1994; Garcia-Bellido & Linde 
1995) the value of the Planck mass before the onset of inflation places an upper bound on the val- 
ues of the inflaton field, above which metric fiuctuations become too large for infiation to occur. 
On the other hand, the Planck mass varies slowly as the infiaton field rolls down the potential, 
and its value remains constant at the end of infiation. Therefore, the value of the Planck mass 
at the end of infiation dictates the formation of structure and quantum physics in thermalized 
regions. By virtue of global stationarity, one can calculate the probability distribution of G and 
the spectrum of density perturbations in those regions in a time-independent manner. In this 
paper we derive the joint probability distributions of the infiaton and BD fields within this sce- 
nario, therefore predicting the likeliest values of the Planck mass at the end of infiation given the 
initial conditions, and the spectrum of density fiuctuations is computed therefrom. 

In order to calculate the spectrum of density fiuctuations in different regions we have mentioned 
that we need to know the volume ratio of these regions. A way to proceed is to choose a foliation 
where the time variable is the synchronous time and calculate the 3-volumes at constant time. 
We construct the ratio 

(I) 

which gives us the relative probability of the infiaton field a taking values in the range ctq, ctq -|- dcr 
with respect to the thermalized regions (cr ^ cr*) at time t. It is expected that this ratio will 
converge to a constant value in the limit t — )• oo. It has been pointed out however (Garcia-Bellido 
et ai. 1994; Linde et ai. 1994,1995) that the calculation of the probabilities is extremely sensitive 
to the choice of time parametrization, and different choices of t yield to inconsistent results in 
the asymptotic limit. 

Recently, a method was proposed (Vilenkin 1995a) whereby one can investigate the ratios of 
physical 3-volumes of different regions of the universe in a way that is very insensitive to the 
choice of time parametrization. The method consists on computing the 3-volumes of hypersur- 
faces (T = const, not at t = const but over all spacetime. Thus, one construct ratios similar to 
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(1), 

<-^' 

where the volumes are integrated for all times. Certainly the individual volumes are infinite, but 
due to global stationarity, the ratio of any two volumes V(cri) and V(cr2) remains finite. The 
procedure to calculate these ratios consists on a regularization proposal that is summarized in 
Section 3. 



From global stationarity, we have that the comoving probability distribution of the infiaton field 
reaches an asymptotic regime, 

V{t^^,a) = Vs{(j), (3) 
which is related to the 3-volumes V via 

VsW) = ^, (4) 

where 

Vt = / d(7 V{(t). (5) 



By calculating the probability distribution of the infiaton field, the regularization method can be 
used to compute the spectrum of density fiuctuations in a way that is insensitive of the choice of 
time variable. 



In this paper, we investigate the application of the regularization procedure of Vilenkin (1995a) 
to models of chaotic infiation in a BD theory of gravity, and we derive the spectrum of density 
fiuctuations. The article is organized as follows: in Section 2 we describe the fundamentals of 
eternal infiation in a gravity theory that is of the BD type for two different potentials; in Section 
3 the regularization procedure for comparing infinite volumes of constant -field hypersurfaces is 
summarized for an infiation-only model (BD field constant); in Section 4 we calculate the prob- 
ability distributions in a generic BD chaotic infiation we apply the tools of Section 3 for the 
derivation of the volume ratios given by these probability distributions; in Section 5 we derive 
the amplitudes of the quantum fiuctuations of the fields measured by a typical observer and we 
use these results to calculate the spectrum of density fiuctuations; finally in Section 6 we discuss 
the results and possible observational predictions. 



2 Eternal inflation in Brans— Dicke gravity 



The BD action in chaotic infiation is given by (e.g. Barrow & Maeda 1991; Liddle & Wands 1992; 
Garcia-Bellido et ai. 1994): 

S = jd\^g [j-^'R - \{d^f - \{daf - Via)] , (6) 
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where R is the spacetime curvature and V{cr) is an arbitrary inflaton potentiaL The Planck mass 
is related to the BD field (f> via 

G-'=MU<P) = -<P\ (7) 
(jj 

where the coupling constant uj satisfies the observational bound uj ^ 500 (Reasenberg et al. 1979; 
Accetta et al. 1990; Casas et al. 1992; Will 1993) and in the limit — )• oo one recovers GR, where 
the value of G is constant throughout. In the context of string theory, the coupling io{(f>) becomes 
a dynamical variable for the so-called dilaton field (f> (see e.g. Bergmann 1968; Nordtvedt 1970; 
Wagoner 1970; Damour & Nordtvedt 1993; Lidsey 1996). The BD field is a particular case of 
dilaton field for which the coupling io with matter is constant, as we shall assume hereafter. 

The Planck boundary or beginning-of-infiation boundary is given by the curve 

V{a) = MU<P), (8) 

that marks the boundary were metric fiuctuations begin to be large, and hence, we do not allow 
the infiaton field to take values for which the potential would be above this boundary. In the 
case of arbitrarily large values of ^, a larger range of values of a is permitted, and thus one may 
encounter the situation where the fields diverge simultaneously while remaining within the curve 
(8). On the other hand, the end of infiation is marked by the condition 

l^' + l^'^Via), (9) 

which, like the Planck boundary (8), is also in principle unbounded. 

The equations of motion in an FRW background become: 

l>' + ^i?)<^ = 0, (10) 

DV = -y'(a), (11) 



^ k S Aid 



where the curvature k = 0, ±1 and the differential operator D is defined 

k^ 



V' = d^ + 3Hdt--. (13) 
a/ 



In the slow-roll approximation, i.e. during infiation, where (f> <^ Hcf) <^ H^(f>^ -\- (f>^ <^ 2V{cr) 
and y"((T), k^ <^ H^^ these equations read 



_ H 

(jJ 

1 



(14) 
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4:LO 



(16) 



and the curvature scalar is given by i? = —12H^. Equations (14)-(16) permit us to write the 
end-ol-inllation boundary (9) in the form 



:i7) 



Thereupon, in the limit of large cu, we obtain that the Planck mass at the end of inflation is given 
by 

W 

(18) 



M|, = 487r 



V 



In our own observable part of the universe, Mp ~ 10^^ GeV, and therefore. 



lO^^GeV, 



(19) 



where cr* denotes the value of the inflaton at the end of inflation. It is easy to show also from 
(14)-(16) that the following general conservation law holds: 



dt 



+ 8 / d(7 



(J 



V'(a) 



0. 



(20) 



Notice also that the evolution of the inflaton field is determined by the BD field through the 
Hubble parameter in the damping term on the RHS of (15). This effect will be greater the 
smaller H is, which from (14) one can see occurs for slow variations of log (p. In the following we 
discuss the dynamics of this scenario for two infiaton potentials. 



2.1 Model I: Power— law potential 

In this case we consider potentials of the type 

Vi<r) = ^ (21) 
Zn 

where n is an integer an the coupling A is typically of order unity. The Planck boundary is then 
given by 

Thus, we require that initially > given an arbitrary ctq. Hence, at the initial time, provided 
the latter condition is satisfied, one may find the situation where both and ctq are very large 
at the onset of infiation but still the infiaton field lies under the Planck boundary. The equations 
of motion in the slow-roll regime become 

(h to 



6 



cr 

(T 



n H (f)^ 



3(^2 2n 

The end-of-inflation boundary is given by the curve 

n 



a 



2n 



2 ^* (3^ - 2) 



(24) 
(25) 

(26) 



where the star-subscript quantities denote the values of the fields at the end of infiation, and 
from (23) (24) again we derive a conserved quantity: 



dt 
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5 H (7 

n 



0, 



(27) 



which tells us that the classical motion of the fields in the slow-roll regime is along a ellipse on the 
(^, cr) plane as is illustrated in Fig. 1. It is apparent from Fig. 1 that, given some initial conditions 
((To,^o), as the infiaton field rolls down the potential during infiation by action of its classical 
drift, the BD field increases in value monotonically and its rate of growth is weakly dependent on 
n. Therefore, the Planck mass at the end of infiation is always greater than its initial value. The 
three figures correspond to potentials n < 2, n = 2 and n > 2. In all cases, the end-of-infiation 
curve is a straight line, whereas the Planck boundary is strongly dependent on the value of n. 
Only in those cases where n > 2 do the Planck boundary and the end-of-infiation boundary 
intersect, and the enclosed finite area is the only region where infiation can take place. This 
effectively sets an upper bound for the values of the fields at the initial times and, ultimately, for 
the value of the Planck mass at the end of infiation. In the other two cases the initial conditions 
are effectively unbounded, provided that they lie in the region above (22). 



Conservation of the quantity <f>'^ + {2/n)a'^ at the initial and end-of-infiation times yields 



2 2 

i2 I ^ 2 ^2 

n n 



I + -(3cu -2) 

n 



Hence, 



(J. 



I + -(3cu - 2) 



(3cu 



I + -(3cu - 2) 



(29) 

If cf>o ^ (To, then the fields at the end-of-infiation boundary are fully determined by the BD field. 
In chaotic infiation, the initial situation is that the infiation field starts out from values close 
to the Planck boundary and rolls down the potential towards smaller values. This implies, as 
pointed out above, that the BD field grows and therefore the Planck mass at the end of infiation 
is larger than its initial value. However, in the case of ^ Co, and in the limit u — )• oo, the 
Planck mass remains practically constant as it is easy to show from (29) (30). We get 



^2 



6cu 



(31) 



7 



^0- 



(32) 



On the other hand, from (18), we get that if a powerlaw potential is a good approximation to 
the inflaton potential in our observable universe, then 



2n 



W GeV, 



(33) 



or equivalently. 



id 



1/2 



10^" GeV, 



(34) 



which requires a large value of (po for most conservative estimates of io. 



2.2 Model II: Double-well potential 

Let us consider the case of a symmetric double- well inflaton potential such as: 

V{a) = \m'{a'-rjy. (35) 
The Planck boundary is then given by the curve 

' ^ ifj^ - ri^\ = <i>\ (36) 



V 47r 

Given cr initially, the Planck boundary is precisely on top of the hill at cr = if 

^^ = ^l = —r,\ (37) 
47r 

with a corresponding Planck mass Mp = (m/2) i]^ . Thus, if cf){t — > 0) > the Planck boundary 
will be located above the maximum of the potential at cr = 0, and if cf){t — > 0) < it will cut 
through the slope of the potential. In the latter case, the allowed initial value of a is greater than 
a threshold value 

1/2 

(38) 



mto 



so that cr will lie under the Planck boundary. The subindex denotes the values of the fields at 
t = 0. The equations of motion in the slow-roll approximation are: 



(7^ fm^ 



2 V ^ ' - (^9) 

cr^ \.iLO , 



2 

\ ' 2 2\2 



3^,X-.r, (40) 



H' = ^i<r'-v')\ (41) 
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and therefore the end-of-inflation boundary is given by 

2 



* \3cu — 2/ ul 

Hence, the inversion of (42) gives us the value of the inflaton at the end of inflation, 



f 



1/2 



2 V3cu - 2 
and the Planck mass at this epoch is 



-1 + 



l + ¥(3--2) 



1/2 



(71 



(42) 



(43) 



(44) 



Thus, in the limit — )• oo, assuming a double- well inflaton potential is a good approximation in 
the observable universe, one gets 

12 21 

fO^^GeV 



(45) 



at the end of inflation and, therefore, a very large value of rj is required such that — ?7| ^ rj. This 
in itself implies that rj ^ f 0^® GeV. Therefore the potential is very flat and inflation terminates 
a long while before a has reached the minimum of the potential. From (39) (40), the conservation 
law (20) takes the form 

d 



At 



+ cr — 1] log cr' 







(46) 



and therefore the motion of the fields is confined to a circle in the (^, u) plane, where u is given 
by 



u 



(T — rj log cr 



(47) 



Let us consider the situation where the infiaton field is initially small, ctq ~ e. The conservation 
law (46) enables us to describe the evolution of the BD field as the infiaton rolls-down the 
potential towards the end-of-infiation boundary. Hence 

4>l - rj' log e'^<Pl + al - rj' log (48) 

where we have neglected a O(e^) term. By imposing the bounds cr^ > and cr^ < we get 



Both bounds on the left and right of (49) diverge in the limit e — )• 0, and thus 



loge 



-oo 



0. 



(49) 



(50) 



Therefore, if the infiaton field starts from cr at the initial time, the BD field will grow 
indefinitely in the slow-roll regime regardless of its initial value (po: 



(51) 
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If however, (1) ctq > (due to the Planck boundary being on the slope of the potential for 
example, i.e. < ^ or (2) cr* is unbounded, then (f> will remain thoroughly finite and thus, 
the resulting Planck mass at the end of infiation will be finite. We can obtain a Planck mass of 
~ 10^^ GeV in our observable universe in the limit — )• oo provided that the ratio ^ is large. 
This will require that either the initial (po is very large, or alternatively, that ctq be close to zero. 
In either case, the possibility of having the Planck boundary across the slope of the potential is 
ruled out because this would imply a small value of ^o, and also a large ctq since, as we have 
concluded above, the potential is very fiat around its maximum and minima due to the large 
value of rj. Thus, a further constraint arises on the parameter m of the potential: 

1/2 

^Pi^*) > (^) ^ > 10'" GeV, (52) 

and therefore m ^ 1, where we have used (c<;/47r)^/^ ^ 10. As is well known (see e.g. Linde 
1990b), we need m ~ 10~^ Mp in order to produce density perturbations over scales of astrophys- 
ical interest Sp/p ~ 10~^, which is consistent with the constraint (52). 



3 Comparing volumes in an eternally inflating universe 



In this section, I summarize the regularization procedure proposed in Vilenkin (1995a) (in an 
infiation-only scenario) for the case of a powerlaw potential (21). In the following section, these 
results will be extrapolated to the BD model, that we will use in order to calculate the spectrum 
of density fiuctuations. At a sufficiently early time, we pick a hypersurface on which a ctq, and 
the radius of the universe is normalized to unity at this epoch, a = 1. Assuming an eternally 
infiating scenario, one is able to choose this normalization at an arbitrary radius of the universe. 
Let us denote the 3-volume of this hypersurface Vq. A number N of observers located on this 
hypersurface will, after a given period of infiation, evolve into the hypersurface cr = cr*, which is 
the value of the field at the end-of-infiation boundary at the slope of the potential. The initial 
volume per observer on the a ^ ctq hypersurface is given by Vi(0) = Vq/N and after a time a*i, 
the z-th observer will be located on the hypersurface cr cr* and its volume will grow by a factor 
of a^-. Therefore, the total 3-volume of the hypersurface cr = cr* is given by 

K = ^E4, (53) 

i 

where the sum is over all the observers that reach cr*. In the limit — )• 00, the sum in (53) 
diverges due to the fact that for any arbitrarily large time there will always be regions undergoing 
infiation, and thus there are a^i terms growing without limit. The regularization procedure that 
we will use consists in setting a cutoff in the sum (53) for a small fraction e of regions with the 
largest values of a^i. Therefore, one can write V* in terms of the parameter e and thus compute 
the ratio of the 3-volumes of any a = const hypersurfaces in the limit e — )• 0. The success of this 
procedure lies in that the regularized volumes of the hypersurfaces of constant a is calculated 
over all spacetime and it is therefore independent on the choice of spacetime foliation. 
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3.1 Volume ratios 



The comoving probability distribution Vcic^ t) of the inflaton field a satisfies the continuity equa- 
tion 

dtV, = -d,J, (54) 



where the probability current 

J = 



OTT"^ 47r 



(55) 



and the parameter a defines a time parametrization t that is related to the proper time r via the 
relation (see e.g. Mijic 1990,1991) 

dt = i7^-"dT, (56) 

where H is the Hubble parameter. The fiux J tells us the fraction of the comoving volume at 
a given hypersurface of constant value of cr, such that |J((T, t)| dt is the fraction of observers on 
such a hypersurface between t and t + dt. The synchronous time parametrization is recovered for 
the case of a = 0, such that a = e* and thus 



V* (regularized) = Vo 
where the cutoff time tc is given by 

J((T*, t) dt 



J(c7*,t) e-"* At 



J((T*, t) dt 



(57) 



(58) 



The cutoff parameter e regulates what fraction of the late-time integrated probability fiux is 
removed from the regularized volume. In the limit e — )• 0, we have V* — > oo. It is easy to show 
that in the regime where diffusion can be neglected, i.e. during the slow-roll regime 



J 



1 H' 
-p 



(59) 



the ratio of the 3-volume hypersurface a = (Ta (where a a is an arbitrary value in the range 
|c*| > ^ 0) with respect to the thermalized volume V* is 



where 



Va 



exp 



1 3 



T<7 0-* J 



_47r / — iild£ 



(60) 



(61) 



denotes the expansion factor during the classical slow-roll from the initial hypersurface a ^ ctq 
and the final cr. The ratio (60) represents the relative probability of finding the infiaton field at 
a value (Ta at an arbitrary spacetime point with respect to the thermalized or end-of-infiation 
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value cr*. By virtue ol global stationarity, this relative probability is time-independent. We can 
rewrite (60) in the form 



r = exp 



127r 



For the case of a double-well potential (35) we get 



row 



(J* 



exp [- 37r(cr^ - cr^)] , 



(62) 



(63) 



and similarly for a powerlaw potential (21) 



rpL = exp 



67r 

n 



(64) 



In the case of the double-well potential, as is well-known, the probability of the field for values 
away from cr* is exponentially suppressed, and rather less so in the case of the powerlaw potential, 
due to the factor preceding the exponential in (63). As shown in Vilenkin (1995a), the quantity r 
is very insensitive to the choice of time parametrization and can be therefore used to good effect to 
compute the spectrum of density fiuctuations, as will be discussed in §5. A mild dependence on the 
choice of parametrization is however present, and this is due to the no-diffusion approximation 
(59). On the other hand, it has been shown recently (Linde & Mezhlumian 1996) that, for 
regions undergoing infiation which are away from the non-diffusion regime, the dependence on 
the choice of time parametrization (as indeed the dependence on the choice of regularization 
procedure) can be rather strong, so the results become ambiguous as we depart from our slow- 
roll approximation. It is expected nonetheless, that in the regime where this approximation is 
valid (away from the Planck boundary and close to the end-of-infiation boundary), the results 
discussed in the remainder of this article will be consistent and insensitive to the choice of a. We 
will discuss the performance of this approximation in the context of BD chaotic infiation in the 
following section. 



4 Volume ratios in Brans— Dicke scenarios 



4.1 Stationary probability distributions 

In this section, we extend the results summarized above for infiation-only scenarios to BD eternal 
infiation. The probability distributions derived apply to generic potentials and we will investigate 
the application to the potentials (21)(35). In BD infiation, the comoving probability 'Pc(^,c, 
of the infiaton and BD fields is given by the conservation equation (see e.g. Starobinsky 1986; 
Linde 1990b) 

dtVc = -daJa - d^J^, (65) 
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where the components of the probabihty current J = (J(j, J^) are defined: 

Ja = -^Hf+' d,{Hf+'Vc) - -^H^-'d^H P„ (66) 

OTT"^ 47r 

= -^Hf+' d^Hf+'V,) - :^H"-'d^H P„ (67) 

and the Hubble parameter H is given by (16). In the no-diffusion [slow-roll] regime that we 
shall adopt hereafter, i.e. on the slopes of the potential, the first term on the RHS of (66) (67) is 
negligible and thus 

J,^-^H"-^d,HV,, (68) 
47r 

J^^-^H"-^d^HV,. (69) 

ZTT 

One can expand the probability distribution Vc{<t>Tcr^t) in terms of eigenfunctions in order to 
transform (65) into an eigenvalue equation. We have 

oo 

P,(<^,a,t) = ^V4<^,^)e-"'"*, (70) 

n=l 

where the eigenvalues 7i are ordered 71 < 72 < 73 < . . .. In the limit t — )• 00 the dominant 
contribution is 

P((^,a,t^ ^) = Vi((^,f7)e-"'i*, (71) 

and therefore (65) for t/ji reads 

- 47r7iVi = [H"-^d,H V-i] + 28^ [H"-^d^H V-i] . (72) 

In the following we calculate the volume ratios for the time parametrization given by a = 0, which 
corresponds to t = log a. As shown in Vilenkin (1995a), the regularization procedure summarized 
in §4 yields results that are insensitive to the choice of a in the slow-roll regime. The eigenvalue 
equation (72) then becomes 

- 47r7i Vi = 5.(Z;Vi) + 2a^(Z;Vi), (73) 
where we have defined for shorthand the quantities 

Z;^ ^ H-'d^H = i^, (74) 

Z:^^=H-^d^H = -r\ (75) 

which only depend on a and cf) respectively. Given the initial conditions (^0, Co), then the solution 
of (73) is 

V'i(^,cr) = 2Co^ (^— ^ exp(T<, + T^), (76) 
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where 

r^a)^-^ r^.(C)dC, (77) 

rM^-^ t (78) 

and Co is a normalization constant. Substituting (74)(75) in these integrals and making use of 
(20), we get 

T,+T^ = TT^l{f -^l), (79) 

regardless of the form of the potential. The term e~'^"'^'^o in (76) can be absorbed in the normal- 
ization coefficient Co as will be assumed hereafter. Hence 

P ^ 2Co exp(7r7i - ^,t). (80) 

The parameters Co and 71 are determined by the normalization of the probability and the form 
of the potential. On the one hand, we have that at the initial time the probability is sharply 
peaked at values of the fields in the neighbourhood of the Planck boundary, thus: 

/ d(jd(f) M(j,(f)) = 1, (81) 

where A/" is a sufficiently broad region neighbouring the Planck boundary. In the case of the 
double-well potential, this is located on the top of the hill in the neighbourhood of cr 0, 
whereas in the case of the powerlaw potential A/" is a region of large values of a and ^, just below 
the Planck boundary. On the other hand, the value of 71 is computed by solving the eigenvalue 
equation for a given potential and magnitude of its coupling constant (see e.g. Linde & Mezh- 
lumian 1993 and Linde 1994b for a similar calculation in the infiation-only scenario). Provided 
that the physical volume of the infiating regions grows as a"^*, from (71) we have that infiation 
will only be eternal if 71 < 3. 

Within the slow-roll approximation, the solution (76) automatically satisfies the boundary con- 
ditions at the end-of-infiation boundary, namely the conditions of conservation of probability 
and probability fiux, since the probability currents employed, (68) (69), are equal at either side of 
the end-of-infiation boundary, the diffusion terms having been neglected. 

In the case of the double-well potential (35), we get 

and for the powerlaw potential (21), 

Vc{(J, <P,t) = Co<p- exp(7r7i - -f^t). (83) 
n 

In both cases, there is a greater probability for regions to continue infiating and for the BD field to 
take large values. Therefore, the Planck boundary is likelier to allow larger values of the infiaton 



cr 



2(7 



exp(7r7i (f)"^ - 7it), 



f82l 
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field in the initial conditions, and the expectation value of the Planck mass is expected to peak 
around a large value. In the case of (82), the probability is higher for the inflaton to remain away 
from cr^ and thus, most regions have values of a on the slope of the potential and therefore 

continue inflating. In the case of (83), the inflaton field is likelier to take small values so that, by 
virtue of the conservation law (27), the BD field will grow as much as possible. 

In order to calculate the probability distribution of the fields along the end-of-infiation boundary 
we substitute (17) in (76) and thus 



= Co [2{3io - 2)] exp(7r7i - j,t), (84) 



which is a monotonically growing function in (f> and its derivative only vanishes at ^ = 0. There- 
fore, the likeliest value of (f> in any scenario is the highest one allowed by the form of the potential. 
The possibility of the so-called "run-away" solutions (i.e. solutions where both fields grow with- 
out limit along the Planck boundary or below it) remains viable in this scenario, unless the 
potential is such that the Planck boundary intersects the end-of-infiation boundary, like in the 
example of powerlaw n > 2 that was discussed in §3. It could be conjectured that if Mp in our 
own neighbourhood has a typical value with respect to other regions in the universe, then the real 
infiaton potential will be such that such an intersection of boundaries ought to take place and 
hence the calculation of Vc would produce a distribution that does not allow runaway solutions. 



4.2 Physical volumes of thermalized and inflating regions 



The total 3-volume V* of thermalized regions in the BD eternal infiation model is determined by 
the two-dimensional probability fiux of the infiaton and dilaton fields across the end-of-infiation 
boundary (43). This probability fiux tells us the fraction of the volume of the universe that has 
undergone thermalization. Along the end-of-infiation boundary, the differential probability fiux 
through a line element dl is given by dl (J - n), where n is a normal vector to the end-of-infiation 
curve (17). Hence we have 

(85) 



where the cutoff time tc is again defined by 



dt e^* / d/ ( J • h] 
Jei 



dt J d/ ( J • h] 



dt / d/ ( J • n] 

Jei 



(86) 



For a given infiaton potential it is easy to show, after a little algebra, that the fiux across the 
end-of-infiation boundary is given by 



EI 



dl (J-h) = — e"'i* / d(j exp 

TT 



87r7i {3co - 2) ( ' 



V_ 

V' 



2(3cu-2) 1 



VV" 



V 



12 



+ 1 



(87) 



This expression may be computed either analytically or numerically depending on the potential, 
but a useful simplification can be introduced by taking the limit — )• oo. The asymptotic form 
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of (87) in this limit then reads 



dl (J ■ n] 



EI 



Co 



exp(7r7i 



lit), 



(88) 



where (pmax is the maximum value the BD field can attain along the end-of-infiation boundary. 
In the case of a powerlaw potential this corresponds to cr* = and in the case of a double-well 
potential = i]^ . Therefore, in this limit (85) becomes 



V* ^ Vo 



Co 



exp 



exp(^7r7i 



(3-7i)^. 



(3-7i) 



and the cutoff time tc is determined by 



— log e. 

7i 



(89) 



(90) 



Similarly, we now proceed to compute the regularized 3-volume V(cr, ^) of the regions with 
arbitrary fields (cr, cf)) still undergoing infiation. In order to do this, we note that the fraction of 
observers located at (cr, cf)) within the lapse of time t,t + dt is given by |(J • /) dt|, where / is the 
tangent vector to the curve that crosses (cr, cf)) given by the conservation law (20): 



Thus 



we nave 



V(a, (^) = Vo 
and therefore from (66)(67)(91) we obtain 



[ V 
+ 8 J da — = const. 

(/•/) 



dt e^* 



(91) 
(92) 



Via, 



Vr 



Co 

47r 



+ 16 



V_ 

V' 



1/2 



exp 



exp(^7r7i( 



(3-7i)^. 



(3-7i) 



Hence, finally the ratio of the volume occupied by regions on a hypersurface (cr, 
to the thermalized regions in the — )• oo limit is then 



V( 



cr. 



27r7i 



+ 16 



V_ 

V' 



1/2 



exp 



7r7i 



maxy' 



(93) 

with respect 
(94) 



or equivalently (via (20)) 
r = 27r7i 



max 



V /y\2ii/2 

8 ' ^--p + 16(^77 



exp(^-87r7i j dcr-^ 



(95) 



Once again we note a tendency of the fields towards large values. For any infiaton potential the 
BD field is most likely to be close to ^max- From the curve (91) this is equivalent to a tendency 
for small values of cr, but this is somewhat dependent on the specific shape of the potential. If 
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(pmax is not bounded, then the ratio r is neghgible for regions with finite values of the BD field. 
Hence, the largest physical volume is occupied by the hypersurface (pmax and if we expect our 
observable universe to be a typical region of the universe, (pmax is likely to be very large but 
finite, (pmax ^ 10^° GeV. The regions where (f> ^^ax are close to the end-of-infiation boundary, 
and therefore the largest fraction of the physical volume is occupied by infiating regions that are 
about to thermalize. For the infiaton field on the other hand, the form of r depends entirely on 
the potential. In the case of a powerlaw potential there is again a tendency towards small values 
at the slope of the potential, in a trade-off between the competing factors {V/V) and e^'^^'^ , 
whereas in the case of the double-well potential r has a sharp peak centred at cr = and it tails 
off at (T^ j]^ . Only in the case of the exponential potential e"'^ is the ratio r insensitive to the 
value of cr. 



4.2.1 Powerlaw potential 



In the case of the powerlaw potential, (87) can be calculated analytically. One obtains 

^lJ-n= ,^ .f° . (1 + ^^) exp(7r7iAa2-7it), (96) 
where 

A = 4(3cu-2), (97) 

and we have retained only the dominant contribution to the integral and used the explicit form 
of the probability current: 

Cq f ^ 2 ^ / 2 



J = (J,,, J^) = — (^^, --crj exp(^7r7i^ -litj- (98) 

Therefore, bearing in mind that the time integration in (85) is dominated by its upper limit, one 
gets 

and similarly 

Ca f 4 expf(3-7i)tJ 
V(a, <^) ^ Vo ^ U + -y) exp(7r7i c^') (100) 

In this case, the conservation law employed tells us that (91) is an ellipse of the form ^^ + (2/n)(T^ = 
const. Thus, the volume ratio r in this case is then 

/ n\~^ f 4 \ 1/2 
r = 7r7i nA [l + -j [cf" + — a^j exp(7r7i(^2 _ T^^f^^al). (101) 

It is easy to see that this expression is consistent with (94) in the — )• 00 limit. The value of ctq 
that dominates the integral (96) is the largest a accessible in the initial conditions. Considering 
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that this value lies under the Planck boundary, it is consistent with the approximation ^^^x ~ ^"^o 
provided that the Planck boundary and the end-of-inflation boundary intersect at ^max, i-e., 
n > 2. The alternative situation would be that of an unbounded ^max and in that case ctq would 
be unbounded too. 



4.2.2 Double-well potential 



In this case, the integral (87) can only be computed numerically, but to illustrate the results in 
an analytical form we use the asymptotic result (88) for the limit — )• oo. Thus V* is given by 
(89), and the volume V(cr, (f>) by (93). Hence, the ratio (94) is then 



27r7i 



+ 



(J 



1] 



2\2 



(7^ 



1/2 



exp 



71-71 1 



maxy' 



(102) 



As was noted in the previous section, this ratio peaks sharply at cr and tails off towards 
(T^ j]^ . Therefore most regions are located in the neighbourhood of the top of the hill in the 
potential, as as we have seen in §2.2, this corresponds to the limit ^max oo- Provided ^max 
remains finite, then the conservation law (46) sets a lower bound for and no regions exist 
where a is exactly zero, but in most regions it is situated in a small neighbourhood around this 
value. 



5 Spectrum of density fluctuations in BD inflation 



In this section, the spectrum of density fiuctuations is calculated following the regularization 
scheme we have used above. We shall first derive the distribution of quantum fiuctuations of 
the fields in BD eternal infiation that is measured by a 'typical' observer at a certain value of 
H. This distribution tells us the expectation value of the quantum fiuctuations, (^8(t) and {Scf), 
that contribute to the classical field over a distance scale . Strictly speaking, the spectrum 
of density fiuctuations is calculated via a Fourier decomposition of the stochastic fiuctuations of 
the fields, such that only those contributions that extend over distance scales contribute 
to the integral. Here we assume that the values of the fields are taken at the time when the 
perturbations enter the horizon and we investigate the fiuctuations for a given value of H for 
simplicity. 

For an ensemble of observers located on a hypersurface (cr^,^^), the fiuctuations of the fields 
follow a Gaussian distribution due to the stochastic nature of infiation, with an additional factor 
that depends on the stage of infiation the fields are in, as discussed below. The spectrum of 
Gaussian fiuctuations is described by the distribution 

dV,{8a,8<i>) = —1— exp(-^^^^) dSadS^, (103) 
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where we take the simple-minded view that the variance of the fluctuations is the same for 
both fields at first order in i7, A ^ (Garcia-Bellido et al. 1994). The fiuctuations around 
the hypersurface (cr^,^^) depend on the physical volume occupied by the hypersurfaces in its 
neighbourhood, as it is likelier that the fields will jump from less probable values to more probable 
ones, and this likelihood is proportional to the ratio of the physical volumes occupied by each 
hypersurface. Thus the probability distribution observed by a typical observer at (cr^, (^a) is given 

dViSa. Si) ~ V(,... + fo.^, + ^») ^ , 1^ I 

V{(Ja, (Pa) 

and therefore the stationary value of this distribution with respect to the variations Scr and 6(f> 
gives us the expectation value of the quantum jumps. These are 



(Sa) 



16 



1 



VV" 



A2 



+ 16 



16 1 



VV 



(105) 



A2 



+ 27r7i 



+ 16 



+ 16 ^7 



(106) 



(27r7iA2 - 1) - A2 



The resulting fiuctuations are not Gaussian. These expressions can be simplified by the following 
approximations. As we have seen in the previous section, the predominant scenario is one where 
^2 ^ 0-2, and the volume occupied by regions with large values of (f> is by far the largest and 
therefore the most typical. Therefore 



VV" 
yi2 



A' 



(107) 



a (p. (108) 

We note that the typical quantum jumps in the BD field get larger as the field grows in the course 
of infiation and, as derived in the volume ratio (94), regions with (f> ^^ax occupy the largest 
volume. It is therefore expected that the typical value of the Planck mass at the end of infiation 
is that which corresponds to (pmax- The amplitude of the fiuctuations of the infiaton field depend 
on the shape of the potential. In the case of the powerlaw potential we have 



(Sa) 



(7 



El 

A2 ' 



(109) 



and for the double-well potential. 



(Sa) 



12^ 



:iio) 



In the former case, the fiuctuations are larger for high a and decrease in the course of infiation, 
and in the latter they are larger at small cr, and similarly they become smaller later as the field 
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rolls down the potential towards ±77. In all cases the quantum jumps of the inflaton field decrease 
for increasing (f> and, as we have seen above, (f> does increase in the course of infiation for generic 
potentials. Therefore, these two factors contribute to the non-Gaussianity of the fiuctuations, 
and the end result is that typical quantum fiuctuations are larger at the earlier stages of infiation, 
and the departure from these typical values is exponentially suppressed. 

In order to compute the spectrum of density fiuctuations given by the expectation values (107)(108) 
we use the standard result (see e.g. Linde 1990b) 

_ 6 ^ (7(5(7 + (j)6(f) 

/o ~ 5 a2 + (^2 ' ^ ^ 

where the perturbations are given in the Einstein frame, ^^j, = <f>'^ g^^^ and they enter the horizon 
during the matter-dominated era. Therefore, within the approximation (107)(108), this yields 

_ffi A _ vv" 



^-P = ^-H^u^ 7^ , (112) 

P 5 H^ + ^H'-' ^ ' 



which in the limit of large ^ is reduced to 

hp 487r (H 



p \H'J M|,i 



:ii3) 



For an arbitrary potential it is easy to calculate the spectrum of fiuctuations given by the super- 
position of all the individual contributions (113) of each mode. It can be safely assumed that, 
given that the BD field is almost constant after the end of infiation and its typical value is ~ ^max, 
then typically cf)^ ~ ^max- On the other hand, as shown in §2, H/ H' at the end of infiation is 
typically one or two orders of magnitude smaller than the Planck constant. This is however not 
satisfactory to make (113) sufficiently small so that, in agreement with observational constraints, 
{Spl p) ^ 10""*, and we require that the infiaton field be much smaller than its typical value at the 
end of infiation. For the case of the powerlaw potential, (7*/n ~ 10^® GeV, and thus one would 
obtain a density contrast {6p/p) ~ 10""* only if predominantly a ~ 10~^(7*. 

In those cases where ^max is unbounded, most of the volume of the universe will be occupied 
by a smooth distribution of matter and the Planck mass will grow without limit. In such scenar- 
ios, regions where the Planck mass at the end of infiation is large but finite are untypical, and 
emerge from quantum jumps across the end-of-infiation boundary that prevent (f> from taking 
arbitrarily large values along the classical trajectory. This behaviour does not arise in the case 
for powerlaw potentials n > 2, for which the end-of-infiation and Planck boundaries cross; these 
potentials permit us to predict a finite value of ^max that tells us the typical value of the Planck 
mass at the end of infiation. 
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6 Conclusions 



In this paper we have investigated the probabihty distributions and spectrum of density fluctua- 
tions in a cosmological scenario of BD theories with eternal inflation. We have especially focused 
on powerlaw and double-well inflaton potentials, giving explicit results in both particular cases. 
In the calculation of the probability distributions, we have studied the solutions of the diffusion 
equation in the slow-roll regime and used these results to compute the ratio of physical volumes 
of hypersurfaces with arbitrary values of the flelds with respect to thermalized regions. In accor- 
dance to the principle of global stationarity, these 3-volume ratios tell us the fraction of physical 
volume occupied by arbitrary values of the flelds, in a time-independent manner, at any stage of 
inflation. 

The physical volumes of the hypersurfaces considered are naturally divergent, and hence we 
have adopted a regularization procedure suggested by Vilenkin (1995a) to compute the ratios of 
any two volumes. Recently it has been shown (Linde & Mezhlumian 1996) that different regu- 
larization procedures can be chosen, that enable us to compute volume ratios in a way that is 
independent of the choice of time parametrization, but which yield mutually incompatible re- 
sults. Previous work on the calculation of the spectrum of density fluctuations in BD inflation 
(see e.g. Garcia-Bellido et al. 1994) shows a strong dependence of the results on the choice of 
time parametrization. On the other hand, the approach undertaken here is dependent on the 
choice of the regularization procedure. This is a fundamental issue to be addressed, but I believe 
that we are moving on the right direction. 

As shown by Winitzki & Vilenkin (1996), the residual mild dependence on the choice of time 
parametrization in the regularization procedure of Vilenkin (1995a) is due to the no-diffusion 
approximation employed in solving the stochastic equations of motion for V. It is shown in the 
same paper that the error brought in by this approximation is of the same order of magnitude 
as the "diffusion terms" neglected (flrst term on the RHS of (55)). In principle, a perturbative 
approach is valid to solve the system up to a given order and restrict the dependence on the 
time parametrization to arbitrary accuracy. At the same time, this perturbative approach would 
miminize the dependence of the results on the regularization procedure, since both effects are 
closely related. Ideally an exact analytical solution of the equations would permit us to calculate 
the spectrum of fluctuations in a way that is fairly independent of the choice of regularization 
employed. 

In §2 we have derived a useful conservation law, (20), that helps us to extract some informa- 
tion of the flnal conflgurations in terms of the initial conditions. In particular, we conclude from 
this conservation law that, for the case of a double-well potential, a scenario where ctq ~ is 
only compatible with (pmax oo- In this case, like in powerlaw potentials n < 2, we flnd that 
'run-away' solutions take place, i.e. the leading contribution to the volume of the universe is 
given by conflgurations where the flelds grow indeflnitely. In the probability distributions we 
have calculated, it is found that in a typical region the BD held has a value that is close to (pmax- 
If (pmax is unbounded, most regions will be totally smooth and cosmic structure will not form. 
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{6p/ p) ~ 0, and thus it will be most unlikely to find regions where is finite. However, such 
scenarios are not entirely ruled out, but would lead us to the conclusion that, if this is a plausible 
scenario in our universe, then we inhabit a rare and highly untypical region. If on the other hand 
our region of the universe is typical, then a realistic infiaton potential ought to yield a finite (pmax- 
If at the initial time ^ Co, ^max is close to (po and therefore its magnitude becomes a particle 
physics problem of initial conditions. 

The conservation law (20) also permits us to write all quantities in the classical evolution in 
terms of one field only. The probability distributions (82) (83) (84) may be written in terms of 
the infiaton field only with the aid of (27) (46), as well as the volume ratios, as given by (95). 
Therefore, one can derive the probability distribution for the density field during infiation, V ^ p. 
For instance, in the case of the powerlaw potential, this is such that the distribution has a sta- 

tionary point at . For a very approximate and conservative estimate of the 

parameters, such as n 2, 71 3, this yields p sa 0.16 A°'^^, which a reasonable typical value to 
seed galaxy formation in our own universe. 

In the study of the typical quantum jumps in §5 one takes into account that both fields evolve 
independently of one another and therefore (20) is not satisfied in that regime. The typical quan- 
tum jumps of the fields during the course of infiation do not follow a Gaussian distribution due 
to biased effect of the volume ratios of the hypersurfaces between which the quantum transition 
takes place. This ratio represents the likelihood of a quantum jump in terms of the fraction of 
physical space occupied by each configuration. The typical quantum jumps take a simple form in 
the limit of large ^, which is, as we have seen from the the probability distributions, the likeliest 
scenario. It is found, as expected, that the amplitude of the fiuctuations decreases in the course 
of infiation and it is also smaller for larger (p. In those regions where (f> grows indefinitely, one can 
to a good approximation consider that the classical evolution is valid throughout and quantum 
jumps are negligible. The spectrum of density fiuctuations obtained in this scenario is expressed 
in terms of the relevant parameters of the potential and ^max, and one is therefore able to test 
models so that, provided ^max remains finite, the values of the parameters are consistent with 
the astrophysical bounds for {6p/p). 
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Figure 1: Field trajectories on the (a, ^) plane for a power-law potential. As discussed in §2.1, 
the classical trajectories are ellipses, their axes determined by the initial conditions and n. The 
dashed line is the Planck boundary and the thick solid line is the end-of-inflation boundary. 
For an initial condition located at A, inflation will occur for as long as the fields roll down the 
potential until they reach B. Top: n < 2; centre: n = 2; bottom: n > 2. 
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Figure 2: Planck boundary (dashed line) and end-of-inflation boundary (thick full line) for the 
double-well potential of §2.2. The allowed values of the fields lie in the region above the Planck 
boundary and infiation only takes place for values below the end-of-infiation boundary. It must 
be noted that there is a small neighbourhood of cr^ rj^ that is never reached as a result of 
infiation but becomes negligibly small in the limit ?y — )• oo. 
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